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(54) Modified DCT signal transforming system. 

(57) An MDCT calculating circuit includes an x01 
calculating circuit for multiplying input signals 
with a forward transforming window and a 
linear forward transforming unit for linear for- 
ward transforming an output signal of the cal- 
culating circuit The linear forward 
transforming unit includes an X02 calculating 
circuit and an calculating circuit for pre- 
processing the output signal of the x 01 calculat- 
ing circuit and an integration and summation 
processing circuit for executing integration and 
summation processing operations on an output 
signal of the pre- processing unit The inte- 
gration and summation processing circuit exec- 
utes an integration and summation operation 
on an N/2 number of input signals from the 
pre-processing unit by grouping a k number of 
input signals as a processing unit and iteratively 
executes the integration and summation pro- 
cessing operations a N/(2*K) number of times 
for outputting a sum total of N/2 number of 
signals. 



O 
O 

m 



LU 



VIVO 
lVH103dS 



0-J0 



Mm 



og! 

e 

I; 



.3i§ 



±ofc 

J3P3 

0-J0 



l_ - 



o 

ro 



















J 















• 








M 








<<* 


5 £ 































M 

X <<E 



O-io 



=2 



IVDIOCnOHOHHD | 



VIVO 3HdWVS 
IVOtDOIONOfcWO 



o 

iZ 



Jouve, 18, rue Saint-Denis, 75001 PARIS 



EP 0 590 790 A2 



This invention relates to a modified MDCT signal transforming device for executing fast calculation of linear 
transformation on digital signals such as speech, audio ar picture signals. 

There has hitherto been known a transform encoding employing a so-called orthogonal transform as a 
type of the high efficiency encoding of encoding chronological sampling data such as audio signals by high- 

5 efficiency bit compression. The transform encoding means encoding after block-based orthogonal transform 
of input signals. Typical of the transform encoding is a discrete cosine transform (DCT). Such transform en- 
coding suffers from a problem of block distortion in which non-continuous block-to-block transitions are per- 
ceived as noises. The conventional practice for coping with such problem has been to have an end of a block 
overlapped with an adjacent block. According to a modified DCT (MDCT), a given block is overlapped with 

w both adjacent blocks each by one-half a block size in such a manner as to preclude superfluous transmission 
of samples of overlapped portions. Thus the MDCT is highly suitable for high efficiency encoding. 

Such MDCT and I MDCT, which is its inverse transformation, are disclosed in "Filter Constraints in Hybrid 
MDCT Containing Plural Block Sizes" by Mochizuki, Yano and Nishiya in SHINGAKU-GIHO, CAS 90-10, 
DSP90-14, pages 55 to 60 and in "Adaptive Transform Encoding with Adaptive Block Lengths Using MDCT" 

is by Uzu, Sugiyama, Iwatare and Nishiya in Extended Abstracts in Spring Meeting of Japan Society of Electronic 
Information Communication 1990, A- 197. The above-mentioned MDCT and IMDCT are explained briefly by 
referring to Fig.1. 

Referring to Fig.1 , a given block of chronological sampling data, such as Jth block, has an overlap with 
the (J-1)th block and the (J+1)th block each by a half block (50% overlap). If the number of samples of the 

20 Jth block is n, N being a natural number, there is an overlap of N/2 samples between the Jth block and the (J 
- 1)th block, while there is similarly an overlap of N/2 samples between the Jth block and the (J +1 )th block. 
Each of N samples of each of these blocks, such as the above-mentioned Jth block, is multiplied by a window 
Wh for forward transform or a pre-processing filter to produce an N number of chronological data 102. 

The characteristics of the pre-processing filter or window Wh for orthogonal transform are selected de- 

25 pending on statistic properties of input signals so that power concentration of the transformed data will be max- 
imum. By processing these N samples of the chronological data 1 02 by linearforward transformation by MDCT, 
an N/2 number of independent spectral data 103, that is one-half the number of input samples of the spectral 
data, are produced on the frequency domain. These N/2 number of the spectral data 103 are processed by 
inverse MDCT for producing an N number of chronological data 104 which are multiplied by a window Wf for 

30 executing inverse transform or a synthetic filter to produce chronological data 1 05. The resulting data are sum- 
med to output results of the forward and backward blocks to restore the original chronological input data. 

Meanwhile, in practical processing operation of MDCT and IMDCT, a method for implementing fast oper- 
ation using fast Fourier transform has been proposed in A research in MDCT System and Fast Operation" by 
Iwatare, Nishiya and Sugiyama, in SHNGAK-GIHO, CAS 90-9, DSP 90-13, pages 49 to 54. The same method 

35 has also been disclosed in JP Patent KOKAI Publication No.4-44099 (1992). The calculating method for MDCT 
and IMDCT disclosed therein is hereinafter explained. 
The definition of MDCT is given by an equation (1). 

40 K ' 1 

Vo <*> =c o£ x c h W cos U (2Jc+l) (2n+JV/2+l) / (2N) ) 



0 < k < N/2-1 * * * ( 1 ) 

45 

In the above equation (1), xo is an MDCT input signal, N a block length, h a window function for forward 
transform, y 0 an MDCT output signal, C 0 a constant, n an integer of from 0 to N-1 , and k an integer of from 0 
to N/2-2. Since MDCT processing on sliced chronological data is executed independently for each of the blocks, 
so the block number J is omitted. Meanwhile, since the value of C 0 is not critical in carrying out the calculation 
of MDCT, the following explanation is made on the assumption that C 0 = 1 . For carrying out the operation, x<> 
is multiplied by a window for forward transform to find x^ as indicated by the equation (2). 

xoi(n) = xc(n)h(n) 0^n^N-1 (2) 
From x^ as found from equation (2), Xo 2 is calculated in accordance with equation (3). 
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30 



35 



45 



-x Q1 (n+3N/4) 0 < n < N/4-1 

x 02 (n) = I x Q1 (n-N/4) N/4 < n < N-1 ■•■(3) 

5 

X03 is calculated in accordance with equation (4). 

x 03 (n) = Xo 2 (2n) - x^N - 1 - 2N) 0 ^ n ^ N/2 - 1 (4) 
Then, X03 is multiplied by a coefficient as shown by the equation (5) to produce a series of complex number 
signals z^. 

10 zo^n) * Xo3(2n)exp( - j27in/N) 0 ^ n - 1 (5) 

This signal series is processed with FFT having a length of N/2 to produce a series of complex number 
signals Zo 2 shown by the equation (6). 

S/2 -1 

is z 02 {k) = Y, z 01 {n) exp{-j2Kkn/ {N/2)} 

0 < k < N/2-1 * • • (6) 

20 

y 0 i is calculated from the produced series of complex number signals, as indicated by the following equa- 
tion (7). 

y 0 i(k) = Re(zo 2 (k)exp( - j2n(k + 1/2)/(2N))) 0 ^ k ^ N/2 - 1 (7) 
The above is coincident with y 0 defined by the equation (1), as demonstrated by the above-cited Research 
25 Concerning MDCT System and Fast calculation. 

On the other hand, definition of I MDCT is given by an equation (8), 



.V/2-1 

x x [n)=C x f(n) £ y. (Jc) cos U (2Jc+l) (2J1+W/2+1) / (2A") ) 



Jr=0 



0 < n < N-1 • • • (8) 



In the above equation (8), is an I MDCT input signal, N a block length, m a block number, f a window 
function for inverse transformation, x t is an I MDCT output signal, C 1 is a constant, n an integer of from 0 to N- 
1 and k an integer of from 0 to N/2-1. Since IMDCT processing on sliced chronological data is executed inde- 
pendently for each of the blocks, the block number J is omitted. Meanwhile, since the value of C 0 is not critical 
40 in carrying out the calculation of MDCT, the following explanation is made on the assumption that C 0 = 1 . For 
carrying out the calculation, y t is re-arrayed in accordance with the following equation (9) to produce y^: 



j y](2k) 0 < k < N/4-i 

y n (k) = I -y 1 (N-1-2k) N/4 < k < N/2-1 

•••(9) 

so y 1t is then multiplied by a coefficient as shown by the equation (1 0) to produce a series of complex number 

signals z u : 

Zn(k) = yn(k)exp( - \2nkJU) 0 ^ k ^ N/2 - 1 (10) 
The complex number signal series is processed by inverse FFT over a length N/2 to produce a series of 
complex number signals z 12 by the equation (11): 

55 
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y/2-i 

5 

0 < n < N/2-1 ■ * • (11) 

The complex number signal series z 12 thus formed is multiplied by a coefficient as shown by the equation 
10 (12) and a real number part is taken out to produce x^: 

Xn(n) = Re(2z 12 (n) exp( - j2rc(n + 1/2)/(2N)) 0 =i n =i N/2 - 1 (12) 
Thisxn is re-arrayed, while changing its sign, as indicated by the equation (13), and is multiplied by inverse 
transforming window. 



15 



20 



*12 



(n) - 



f (nJXudi+Ny^) 0 £ n £ N/4-1 

-f (n)xll(3N/4-l-n) N/4 S n i 3N/4-1 . ( 13 ) 

-f (n)x u (n-3N/4) 3N/4 <; n <: N/2-1 



The above is coincident with y 0 defined by the equation (8), as demonstrated by the above-cited Research 
Concerning MDCT System and Fast calculation. 

However, since the fast processing of MDCT for finding the spectrum for each domain having a length N 
and I MDCT which is its inverse operation are carried out with the present method using a complex number 

25 FFT having a length N/2, the N/4 log(N/2) times of multiplication of complex numbers, the N/4 log(N/2) times 
of addition of complex numbers and a work area for storing N/2 number of complex numbers are required if 
the algorithm as described above all in pp. 196 to 198 of "Fast Fourier Transform", by E. ORAN BRINGHAM, 
translated by Miyagawa and Imai, Theory of Base 2 FFT Algorithm, is used, meaning that the method cannot 
be said to be sufficiently effective if it is attempted to carry out fast processing by a small-scale processing 

30 unit. 

Various aspects of the invention are set out in claims 1 and 8 appended hereto. Preferred features are set 
out in the dependent claims. 

An embodiment of the present invention can provide a signal processing method by modified DCT in which 
the operation of the modified DCT and I MDCT, its inverse transform, may be realized by integration and sum- 
35 mation based on a K number of terms as a processing unit with the use of a small-sized processing unit. 

An embodiment of the invention provides a modified DCT signal transforming device for processing an N 
number of sample input signals, comprising a forward transforming window unit for multiplying input signals 
with a forward transforming window and a linear forward transforming unit for linear forward transforming out- 
put signals of said forward transforming window unit, characterized in that the linear forward transforming unit 
40 comprises a pre-processing section for pre-processing the output signals of the forward transforming window 
unit and an integration and summation processing section for executing integrating and summation processing 
operations on output signals of the pre-processing section, the integration and summation processing section 
executing integration and summation processing operations on an N/2 number of input signals supplied from 
the pre-processing section, the N/2 number of the input signals being previously grouped into sets each con- 
45 sisting of a K number of terms, said integration and summation operations being iteratively executed an N/(2*K) 
number of times for outputting a sum total of N/2 number of signals. 

It is noted that K in the convolution unit is a power of 2 and is not less than 4. 

The pre-processing unit outputs, with an N number of samples of input signals of from 0 to (N-1), obtained 
on multiplying the forward transforming windows in the forward transforming window section, the n + 3N/4th 
so input signal inverted in sign as the nth intermediate signal for n of from 0 to N/4-1, and an n - N/4th input signal 
as the nth intermediate signal for n of from N/4 to N-1. With an N number of samples of intermediate signals 
of from 0 to (N-1), the pre-processing unit outputs the 2n-th intermediate signal less the N-1-2n-th intermediate 
signal as an nth output signal. 

The integration and summation processing unit decreases the number of times of the processing opera- 
55 tions by grouping an N/2 number of input signals supplied from the pre-processing unit into sets each consisting 
of a K number of terms, based on periodicity, for outputting the N/2 number of the signals. 

The integration and summation processing unit performs, on an N/2 number of input signals supplied from 
the pre-processing unit, a series of processing operations for m of from 0 to N/(2K)-1 until the number of output 
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signals becomes equal to N/2. The series of the processing operations consists in providing, as a jth first in- 
termediate signal, a sum for Is of from 0 to N/2-1 of a product of the (Kl+j)th input signal multiplied by 
cos(rt(2m+1)(4K+1)/(2N), for j of from 0 to K-1 , providing, as a jth first intermediate signal, a sum for Is of from 
0 to N/2-1 of a product of the (K(1-1)+j)th input signal multiplied by sin(n(2m+1)(4KI+1)/(2N), for j of from K to 
5 (2K-1), providing, as a jth second intermediate signal, a sum of the jth first intermediate signal multiplied by - 
cos(2n(2m+1)(j/N) and the (j-K)th first intermediate signal multiplied by -sin(2n(2m+1)(j/N), for j of from K to 
K-1, providing, as a jth second intermediate signal, a sum of the (j-K)th first intermediate signal multiplied by 
-sin(2it(2m+1)(j-K)/N) and the jth first intermediate signal multiplied by -cos(2n(2m+1)(j-K)/N), for j of from K 
to 2K-1, providing, as a k-th third intermediate signal, a sum of the jth second intermediate signal multiplied 
10 by cos(2rcJk(4j+1)/(4K)), for j of from 0 to K-1, providing, as the k-th third intermediate signal, a sum of the jth 
second intermediate signal multiplied by sin(2uJ<(4(j-K)+1)/(4K)), for k of from K to (2K-1), and providing, as a 
(m+kN/(2K))th output signal, a sum of the k-th third intermediate signal and a (m+kN/(2K))th third intermediate 
signal, for k of from 0 to K-1. 

The integration and summation processing unit performs a processing operation of computing an output 
15 signal of the integration and summation processing unit from the second intermediate signal by a processing 
volume reduced by one half by grouping the second intermediate signals by each two terms by taking advan- 
tage of periodicity so that the number of output signals becomes equal to N/2. 

The integration and summation processing unit performs a processing operation of computing an output 
signal of the integration and summation processing unit from the second intermediate signal by a processing 
20 volume reduced by one quarter by grouping the second intermediate signals each by four terms by taking ad- 
vantage of periodicity so that the number of output signals becomes equal to N/2. 

A modified DCT signal processing device for processing an N/2 number of independent input signals, com- 
prising a linear inverse transforming unit for linear inverse transform of input signals, and an inverse transform- 
ing window unit for multiplying an output signal of the linear inverse transform unit by an inverse transforming 
25 window, wherein the linear inverse transforming unit comprises a pre-processing unit for pre-processing an 
input signal , an integration and summation processing unit for executing an integration and summation proc- 
essing operations on an output signal of said pre-processing unit, and a post- processing unit for post- 
processing an output signal of the integration and summation processing unit, the integration and summation 
processing unit executing the integration and summation operations on the N/2 number of input signals sup- 
30 plied from the pre-processing unit by grouping the input signal by a K number of terms, and executing the op- 
erations iteratively by N/(2*K) number of times for outputting a sum total of an N/2 number of signals. 

In the integration and summation processing section, K is a power of 2 which is not less than 4. 

The pre-processing unit outputs a signal equivalent to a k-th output signal which is the 2k-th input signal 
for values of k of from 0 to N/4-1 and a signal equivalent to a k-th output signal which is the <N-1-2k)th input 
35 signal inverted in sign for values of k of from N/4 to N/2-1 . 

The integration and summation processing unit decreases the number of times of processing operations 
for integration and summation processing of N/2 number of input signals supplied from the pre-processing unit 
based on periodicity for outputting said N/2 number of signals. 

The integration and summation processing unit performs, on an N/2 number of input signals supplied from 
40 the pre-processing unit, a series of processing operations for m of from 0 to N/(2K)-1 until the number of output 
signals becomes equal to N/2. The series of the processing operations consists in providing, as a jth first in- 
termediate signal, a sum for Is of from 0 to N/2-1 of a product of the (Kl+j)th input signal multiplied by 
cos(7i(2m+1)(4K+1)/(2N), for j of from 0 to K-1, providing, as a jth first intermediate signal, a sum for Is of from 
0 to N/2-1 of a product of the (K(l-1)+j)th input signal multiplied by sin(rc(2m+1)(4KI+1)/(2N), for j of from K to 
45 (2K-1), providing, as a jth second intermediate signal, a sum of the jth first intermediate signal multiplied by - 
cos(2n(2m+1)(j/N) and the G-K)th first intermediate signal multiplied by -sin(2n(2m+1)(j/N), for j of from K to 
K-1, providing, as a jth second intermediate signal, a sum of the 0-K)th first intermediate signal multiplied by 
-sin(2jc(2m+1)(j-K)/N) and the jth first intermediate signal multiplied by -cos(2n(2m+1)(j-K)/N), for j of from K 
to 2K-1, providing, as a k-th third intermediate signal, a sum of the jth second intermediate signal multiplied 
so by cos(2nJk(4j+1)/(4K)), for j of from 0 to K-1 , providing, as the k-th third intermediate signal, a sum of the jth 
second intermediate signal multiplied by sin(2itk(4(j-K)+1)/(4K)), for k of from K to (2K-1), and providing, as a 
(m+kN/(2K))th output signal, a sum of the k-th third intermediate signal and a (m+kN/(2K))th third intermediate 
signal, for k of from 0 to K-1. 

The integration and summation processing unit performs a processing operation of computing an output 
55 signal of the integration and summation processing unit from the second intermediate signal by a processing 
volume reduced by one half by grouping the second intermediate signals by each two terms by taking advan- 
tage of periodicity so that the number of output signals becomes equal to N/2. 

The integration and summation processing unit performs a processing operation of computing an output 
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signal of the integration and summation processing unit from the second intermediate signal by a processing 
volume reduced by one quarter by grouping the second intermediate signals by each four terms by taking ad- 
vantage of periodicity so that the number of output signals becomes equal to N/2. 

The pre-processing unit outputs a signal equivalent to an (n+N/4)th input signal for the value of n of from 
5 0 to N/4-1 , a signal equivalent to a 3N/4-1-nth Input signal reversed in sign for values of n of from N/4 to 3N/4- 
1 and a signal equivalent to an n-3N/4th input signal reversed in sign for values of n of from 3N/4-1 to N/4-1. 

In an embodiment of the invention, the integration and summation processing unit in the linear forward 
transforming unit or in the linear inverse transforming unit performs an integration and summation processing 
operations on an N/2 number of input signals outputted from the pre-processing unit by grouping the input sig- 
10 nals by a K number of terms and performing the operations iteratively by an N/(2*K) number of times for pro- 
ducing a sum total of N/2 signals, for simplifying the constitution. 

The invention will be described hereinafter, by way of example only, with reference to the accompanying 
drawings, in which: 

Fig.1 is a diagrammatic view showing the operating sequence of modified DCT or MDCT operations and 
15 an inverse MDCT or IMDCT operations. 

Fig. 2 is a block circuit diagram showing an arrangement for realizing an MDCT signal transforming method 
according to an illustrated embodiment. 

Fig. 3 is a block circuit diagram showing an arrangement for carrying out the IMDCT signal transforming 
method according to an illustrated embodiment 
20 Fig .4 is a flow chart showing the process steps of a basic embodiment of an MDCT signal transforming 

method according to an illustrated embodiment. 

Fig.5 is a flow chart showing the process steps of a basic embodiment of an IMDCT signal transforming 
method according to an illustrated embodiment. 

Fig. 6 is a block circuit diagram showing an example of a circuit of a high efficiency encoding device for 
25 carrying out the MDCT signal transforming method according to the illustrated embodiment 

Fig. 7 is a block circuit diagram showing an example of a circuit of a high efficiency decoding device for 
carrying out the IMDCT signal transforming method according to the illustrated embodiment. 

Referring to the drawings, preferred embodiments of the signal transforming device by the modified DCT 
is explained in detail. 

30 Specifically, the signal transforming device by MDCT of the illustrated embodiment (MDCT processing cir- 
cuit 30) has an arrangement as shown in Fig. 2, while the signal transforming device by IMDCT (IMDCT proc- 
essing circuit 50) has an arrangement as shown in Fig. 3. Before proceeding to description of the signal trans- 
forming devices by MDCT and IMDCT of the embodiments shown in Figs.2 and 3, the basic principle of the 
signal processing method (signal transforming method) applied to these devices is explained. 

35 With the signal processing method applied to the device of the present embodiment, a coefficient table 

shown by the following equations (14) to (16) is provided for an integer m which is not less than 0 and not more 
than (N/(2K)-1): 

a mil = cos(7i(2n + 1)(4K + 1)/(2N)) 0 ^ 1 ^ N/2 - 1 (14) 
b Wil = sin(7i(2m + 1)(4K1 + 1)/(2N)) 
40 c mt1 = cos (2Ti(2m + 1 )j/(N) 0 ^ j ^ K - 1 (15) 

= sin(2K(2m + 1)j/N) 
e kJ = cos (2nk(4j + 1 )/(4K)) 0 ^ j ^ K/4 - 1 (1 6) 
f kJ = sin(2nk(4j + 1)/(4K)) 
For integers m which is not less than 0 and no more than (N/(2K)-1), the following calculation is made re- 
45 peatedly. 

First, using the coefficients of the equation (14), a K number of signals A01, as given by the equation (17), 
and a K number of signals Ao 2 , as given by the equation (18), are calculated from x 03 of the equation (4). 

Nf (2f0 -1 
1=0 

0 < j < K-1 • ■ • ( 17) 

55 



6 



EP0 590 790 A2 



0<j<K-1 - ■ ■ ( IB) 

Then, from the coefficient of the equation (15) and from these signals Aql Aq 2 , two signal trains or rows 
10 B 0 i, B 02 , each being of a length K, are calculated as shown by the equations (19) and (20). 

B m ,oi(j) 55 c mJ V01O) - dmA^G) 0 S J s K - 1 (19) 
B m , 02 G) = • dmj Kwti) • Cbg A m , 02 (j) O^j^K- 1 (20) 
Using the coefficients of the equation (16) and B 0 i, B 0 2, a K number of signals Xo are obtained, as shown 
by the equation (21). 

15 



20 



25 



30 



35 



40 



45 



50 



55 



X, (sn+kN/ (2J0 ) =£ ej ,,S m , 0l (j) f kij B m _ 02 (j) 



0 1 m < N/(2K)-1 

0 < k < K-1 • • • (21 ) 

This is coincident with y 0 as y 0 defined by the equation (1), as may be seen from the following: 
First, from the relationship of the equations (1), (2), (3) and (4), the equation (22) holds: 



.--1 

y c (k) =£* 01 (i3)cos(2*;(2i:+l) Un+N+2 ) / ( BN) ) 

j:-0 



5.V/4-1 

£ *oi(n-.V/4)cos(2*(2Jr+l) (4 (n-A r /4) +A>2) / (8A 7 ) ) 
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= £ x cl (n-^/4)cos(2z(2Jc+l) (2/2+1) / (4/7) ) 



5.V/4-1 

+ Y, x Q1 (n-N/4)cos{2%(2k+l) (2/2+1)/ (4/7)) 



*v-i 

= 53 x 0I (n-N/4)cos(2K(2ic^l) (2/2+D/UW) ) 



*V/4-l 

+ x Q1 (n+3^/4) cos (27v (2Jc+l) (2 (n+k) + 1) / (4/7) ) 



iV-1 

= 53 x 01 (n-iV/4)cos (2ti (2Jc+l) (2/2+1)7(4/7) ) 



- 53 Jf 0; (i3+3N/4)cos(2T: (2ic+l} (2/2+1) / (427) ) 



.v-i 

=53 x 02 < n > cos (21C (2ic+l) (2n+l) / (4/7) ) 



N/2-1 

= 53 x 32 (2/2)cos(2ts (2Jc+l) (4/2+1) / (4A 7 ) ) 



*v/2-l 

+ 53 x c2 (2n+l)cos (2* (2Jc+l) (2 (2/2+1) +1) / (4/7) ) 



.V/2-1 

= 53 x Q2 {2n)cos{2r. {2k+l) (4/2+1)7(4/7) ) 

n=0 
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AT/2-1 

r Y, x c2 (W-l-2n)cos(2*(2Jc+l) (2/ML-4ji) / (4W) ) 

2 = 0 



= £ x 02 (2n)cos(2iu(2J:+l) (4n+l)/(4W)J 

n=0 



*V/2-l 

n=0 



>r/2 -1 

5 £ U oa (2n)-x oa (W-l-2n))cos(2ic(2Jc+l) (4n+l)/(4W) ) 

n-0 



*V/2-l 



- ■ ■ (22) 

Modifying this, the equation (23) is derived. 

y 0 0n+JcN/(2f0) = £ x <* ( 1] cos (7r (rn+2 {kN/ {2K) ] +1) (4 J + 1) ' (2AT) 5 

1-0 



X 03 (J) COS (it (2/77+1) (4i+l)/(2iV)+27iJc(42+l)/4in 

i=0 



= E E -^03 CiTi+j) cos {n (2ia+l) (4U2+j)+l)/(2,V) 



+ 2uic(4 (ffi+j) +1) / (4i0 ) 
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r-i K/C2JO-1 

£ x 03 (iCj+i)cos(7i (2/n+l) (4JCJ+1)/(2JW 

j*0 1*0 

+2n (2J7I+1) j/^+2icJc(4j+l) / (4i0 ) 
= V cos (2n {2;n+l) j7tf+2iric(4j+l) (4tf) ) 

x £ (*2 + J) cos (2m+l) (4JC2+1) / (2A0 ) 

1=0 

sin(27c (2m+l) j/N+2nk(4j + l) / UK) ) 
x £ x C3 (Kl+j) sin (i; (2jh+1) (4*2+1) / (2N) ) 

1*0 

0 < k < K-1 •••(23) 
Using the equations (14), (15), (17) and (18), the equation (24) is derived. 

y Q {m+kN/ (2/0 > =£ cos (2- (2/n+l) j/N+2nk{ (4 j+1) / UK) )-!a, #c1 (j) 

K-l 

sin(2- (2/71+1) j/N+2xk( (4j+l) / (4/0 M m , 02 (j) 
=J2 cos (275 (2/n+l) j'AV) cos (27iic (4j+l) / (4 JO M„ #C1 (j) 
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sin(2TC(2m+l) j/M sin(27iic(4 / (4 JO }^ cl (j) 



sin (2ti (2OT+1) j/iV) cos (2k*(4 / (4JO )^ c2 (j) 



cos (25i(2m+l) j'Av) sin(27iic(4j+l) / (4i0 >An,c2 0') 
x-i 

= E c .-n,j cos (2icic(4j+l) / (4K) U a , 01 (J) 

.T-l 

£ d^sinUicJcUi+l) / (4JQ )A at01 ij) 

"E d *.J cos (2wic(4j+l) / (4 JO ) A a>02 (J) 
r-i 

~E c aiJ sin(2ici:(4j+l)/(4J0 M fl>32 (j) 
cos(2icJc(4j+l)/(4J0 ) (c^, 01 (j) -d^ 5(02 (j) ) 



sin (2iE*(4j+l) / UK) ) ( -cf fli ^ 01 (j) ~c rttj A mt02 (j) ) 



cos (2ic*( (4 / (4J0 ) B 9t0l (j) 



sin(2icJc(4j+l) / (4Jf) ) B fli02 {j) 

0 < k < K-1 - - . (24) 

11 
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From the equations (19) and (20), the equation (24) is changed to the equation (25). 



10 



20 



25 



30 



40 



45 



y 0 Un+kN/ (2 JO ) =£ e kij B mi , x (j) + £ f ki ,B a>02 ( j) ) 

J*0 j'O 



0 < m <N/(2K)-1, 0 < k < K-1 — (25) 



This coincides with the equation (21) for the method embodying the present invention. Itfollows from this 
that the MDCT processing may be achieved with the method embodying the present invention. 

If the periodicity of the equation (16) is utilized, the amount of convolution processing of the equation (21) 
may be reduced to one half or to a quarter. 
is That is, the coefficient of the equation (16) is modified to the form shown in Equations (26) and (27). 

e 2k+p,j+K/2 = cos (2TT(2k+p)(4(j + K/2) + 1)/(4K)) 
= cos (2it(2k + p)(4j+l )/(4K) + (2k+p)n) 



= ("^Skipj P = O f 1 (26) 

f 2k+P.j+K/2 = s1n (2TT(2k+p)<4( j + K/2)+l )/<4K) ) 



= sin (2Tr(2k+p) (4j + l )/(4K) + (2k+p)Tr) 

= ( - 1)Pf 2k + p,j P = 0,1 ... (27) 



Therefore, as shown by equation (28), by setting, as shown by the equation (28), 

E m ,o0) = B m ,oi(j) + B m , 01 (j + K/2) 
Em.i(D = B m , 01 (j) - B m ,oi(j + «/2) 
35 F m , 0 (j) = B mfO2 0) + B mtQ2 0 + K/2) 

F m ,i(j) = B m(02 (j) - B m , 02 (j + K/2) 
0 ^ j ^ K/2 - 1 (28) 

Xq is expressed as shown by the equation (29), so that the amount of integration and summation processing 
operations is halved. 



50 (31): 



X/2-1 K/2-1 

X a {m+(2k+p)N/(21Q) = Z e u.pApW) + E W^.pO') 

P = 0, 1 0 < m < N/(2K)-1, 0 < k < K/2-1 --(29) 
Besides, since the coefficient of the equation (16) may also be expressed as shown by the equations (30), 



55 



12 



EP 0 590 790 A2 



e 4k+p,j+qk/4 = cos(2Tr(4k+p) (4( j+qK/4)+1 )/(4K) ) 
= cos(2TC(4k + p)(4j+1 )/(4K) + pqTT/2.)! 
= cos( pqir/2 )e 4k+P(j -si n( pqA/2 )f 4k+Pf j 

p = 0, 1 ,2,3, q = 0,1,2,3 ... (30) 

10 

f 4k+p, j+qk/4 = sin(2K(4k+p)(4( j+qK/4)+1 )/(4K)) 
= sin(2TT(4k+p) (4j + l )/(4K) + pqTr/2) 

15 = Sin ( P^ Tr / 2 ) e 4k+pJ~ sin (P^ X / 2 ) f 4k+pJ 

p = 0,1 , 2 , 3 , q = 0,1,2,3 '...(31) 

20 the equation (32) holds. 

25 
30 
35 
40 
45 
50 
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10 



15 



20 



30 



35 



40 



45 



/r/4-1 3 
j=o <r»o 



*/4-i 3 



-sin (pgrc/2 ) f 4Jt ^) 5 a . 01 {j+qK/ 4 ) 



*/4-i 3 

V) £ (sin(pgn/2)e 4jetptJ 



J*0 cpo 



25 + co s ( pg* / 2 ) f^ pj ) B Bt 02 (j + qK/4) 



K/4-1 3 



■sin(pgn/2)S a 02 (j+gK/4) 



*/4-l 3 

~ £ f «*♦*.*£ (sin(pgn/2)B a , 01 (j+gJir/4) 

J=0 <no 



+cos(pgrrc/2)3 J! , 02 (j'+gi?/4) ) 

0 < m < N/(2K)-1, 0 < k < K/4-1 , p = 0, 1, 2, 3 

• • • (32 



Since one of sine and cosine values in the equation (31) is S1 with the other value being 0, by setting so 
that, as shown in equation (33) and equation (34) 
so G m , 0 (j) = B m , 01 (j) ♦ B m , 0 ,0 + KM) + B m>ol 0 + 2K/4) + M m . o ,0 + 3K/4) 

G mi1 (j) = B m ,oiO) + B m , O2 0 + K/4) - B m , 0 ,(j + 2K/4) - M„, 01 (j + 3K/4) 
G m , 2 (j) = B m , 0 i(i) - B m , 01 (j + K/4) - B„,, 01 (j + 2k/4) + M^O + 3K/4) 
G m , 3 0) = B m , 01 a) - B m . O2 0 + K/4) - B,„ 01 (j + 2K/4) + M^O + 3K/4) 
G m .o0 + K/4) = B m ,o2(j) + B m .o2(j + K/4) + B mf02 0 ♦ 2K/4) + M m , 02 (j + 3K/4) 
« G m>1 (j + K/4) = B^ojO) - B m>a2 (j + K/4) - B m , 02 (j + 2K/4) + M m , 01 (j + 3K/4) 

G m>2 (j + K/4) = B^raO) - B m>02 (j + K/4) + B m>02 (j + 2K/4) - M m , 02 (j + 3K/4) 
G m .,0 + K/4) = B m , 02 (j) + B m , 01 (j + K/4) - B mfi2 Q + 2K/4) - M m . 02 (j + 3K/4) 0 S j S K/4 - 1 (33) 
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j COS(2TTk(4j+1 )/(4K)) 0 < j < K/4-1 

g kj = I sin(2TTk(4( j-K/4) + 1 )/(4K)) K/4 <j < K/2-1 --(34) 

5 

the equation (35) 

K/2-1 

10 Jf 0 (m+(4Jr+p)^V/ (2 JO ) = £ ^J^A,*^ 

0 < k < K/4-1, p = 0, 1, 2, 3 ---(35) 

f5 

holds, so that the volume of the integration and summation operation of the equation (21) may be halved. 

Besides, in carrying out the operation of IMDCT with the signal processing method of the present inven- 
tion, coefficient tables of the equations (14), (15) and (16), which are the same as those for MDCT, are used 
for integers m which are not less than 0 and not more than (N/(2K)-1). The following operations are repeated 
20 for integers m which are not less than 0 and not more than (N/(2K)-1). 

First, using the coefficients of the equation (14), a K number of signals P 0 i as given by the equation (36) 

^,oi (J) = £ a^^iKl+j) 

25 i=0 

0 < j < K-1 • • - (36) 
30 and a K number of signals P 0 2 as given by the equation (37) 



N/ {210-1 

35 1=0 

0 < j < K-1 . . • (37) 

40 are calculated from y^ of the equation (9). 

Then, from the coefficients of equation (15) and from these signals P 0 i, Pq2. two signal trains Q01, Q02. 
each having a length K, are calculated, as shown by the following equations (38), (39): 
Qm.oi(j) = WV01O) - d mJ A m(02 (j) 0sij^K-1 (38) 

Qm.020) = " d m jP m ,oiO) " C m,/W(j) 0 ^ j ^ K - 1 (39) 

45 Then, using the coefficients of equation (16) and B 0 i, B02, a K number of signals W 0 are obtained, as shown 

by the equation (40): 

W 0 (m+kL/ (2K) ) e ktJ Q m01 ( j) f k .fi mM U) 

0 < m < N/(2K)-1, 0 < k < K-1 --(40) 

55 

These signals W 0 , obtained by the equation (40), are processed with a sign change, re-arraying and win- 
dowing for inverse transform, as shown by the equation (41), for calculating an N number of 
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f (n)W 0 (n+N/4) 0 < n < N/4-1 

W^n) = [ -f (n)W Q (3N/4-1-n) N/4 < n < 3N/4-1 
-f (n)W Q (n-3N/4) 3N/4 < n < N/1 

• - -(41) 

This may be found to be coincident with x 1 defined by equation (8) as shown by the following. 
First, an equation (42) holds from the equations (8) and (9). 

.V/2-1 

Xiin-N/i) =C x f (n) £ y, (*) cos (2ic (Jr+l/2) (n+1/2) /N) 

Jt-0 



.V/4-1 

= qjf(n) Y, y(2i)cos(2ic(2J+l/2) (n+l/2)/N) 

2=0 



- Y, y(2i+l)cos(2TC(2J+l+l/2) (jj+1/2) /A 7 ) 

1-0 



= c i f (^) £ y(2i)cos(27c(2J+l/2) (n+l/2)/N) 

J=0 



.V/2-1 

- £ y(2i+l)cos(27:(//-l-2J) /N) 

1 =*V/4 
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= qf(/2) £ y(22)cos(2u(22+l/2) {n+l/2)/m 

1-0 

y/2-1 

" Y, y^-l-22)cos(2* (tf-22-1/2) (n+l/2)/N) 
=qf(j2) £ y(22)cos(27c(22+l/2) (n+l/2)/M 

1=0 

AV2-1 

- £ y(W-l-22)cos(2ic(-2i-l/2) (n+l/2)/W 

.V/4-1 

=qf(n) ^ y(2J) cos (2tc (21+1/2) (j2+i/2)/A-) 

A*/2 -1 

- £ y(iV-l-22)cos(27:(2i+i/2) (n+1/2) /N) 

*V/4-l 

=qf{n) 53 y 11 (J)cos(27t(2i+i/2) (a+i/2) /N) 

1 3 0 

A72-1 

" £ yn(^)cos(27u(2J+l/2) (n+l/2)/A7) 

AV2-1 

^/(js) £ y lx (i) cos (2-n: (2 J + l/2) (n+1/2) /N) 

i*0 



If the equation (43) 

://2-i 

u(n)= £ y xl (2) cos (2* (22+1/2) (n+l/2)/«) 
i-o 
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is set, an equation (44) 



•r-i 

u(n+kN/{2K) ) cos (2itJc(4j + l) / (4 JO ) O at01 ( j) 



sin {2nkUj+l) I (4 JO ) £? flf 02 <J> 



. . . (44) 

holds from the equations (14), (15), (36), (37), (38) and (39), in the same manner as the equation (24). 
Since the equations (45) and (46) 



y/2-1 

u(n+N) = £ y,,{2)cos(27i (22+1/2) [a+N+1/2) /m 

i-0 



= £ y 1 i(J)cos(2i5(2i + l/2) (n+l/2)/tf+(42+l)w) 

i=0 



- X) yn(J)cos(2ic(2i+l/2) (ji+W+1/2) / A 7 ) 



= -u(n) • . • (45) 



.v/a-i 

u<tf/2-l-J3) = J] y^UJcos (2ic (22+1/2) (N/2 -1-J3+1/2 ) /W) 



AT/ 2-1 

= X) yn(i)cos(2ic(2i+l/2) {N/2-n-l/2) /N) 

i«0 



= £ y,, (2) cos (2tc (22+1/2) (-A~/2+n+l/2) /JV) 



1=0 



= u(-N/2+n) 

= -u(N/2+n) . . . (46) 



hold from the equation (43), the equation (47) 
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u{n+kN/ (2 JO ) e^£> m , ci ( j) f*. (j) 

5 

0 < m < N/(2K)-1 t 0 < k < K-1 .--(47) 

holds on comparison with the equation (41), so that it is seen that th IMDCT calculation may be realized with 
10 the method of the illustrated embodiment of the present invention. 

It is noted that, as in the case of MDCT, the volume of the integration and summation operations of the 
equation (40) may be reduced to one half or to one quarter by taking advantage of the periodicity of the equation 
(16). 

First, if equations (48) 
15 S m , 0 (j) = Q m ,oi(i) + Qm,oi(j + K/2) 

S m ,i(j) = Q m .oi(j)-Gm,oiG + K/2) 

T mJ ofl) = Q m ,02(j) + Qm,O20 + K/2) 

T m .i(j) = Q mt02 (j) - Q m ,o2(j + K/2) 
0 ^ j =S K/2 - 1 (48) 

20 are set, an equation (49) 



A72-1 fT/2-1 



25 



W 0 (m+{2k+p)N/{2K))=Y, W*p(J 1+ I] WiW^ 

0 < m < N/(2K)-1, 0 < k < K/2-1, p = 0, 1 --(49) 

30 holds, if the equations (26), (27) are used, so that the volume of integration and summation operations of the 
equation (40) may be reduced by one half. 

On the other hand, if the equations (30), (31) are used, an equation (50) 

35 -V/4-1 3 

u{n+{4Jc+p)iv7(2J0)= J eik*p.iY, (cos (P<?*/2} Q MtQl ( 



40 -sin(pgn/2) 02 (j+qK/4) ) 

.V/4-1 3 

"E f ^P,jY, (sirx(pgn/2)j? a#01 (j + gJC/4) 



45 



50 



+cos (pqn/2)Q sn 02 ) 

O < m < N/(2K)-1, 0 < k < K/4-1, p = 0, 1, 2, 3 

• * - (50) 



55 holds. Since one of sine and cosine values in the equation (50) is ^1, with the other value being 0, by setting 
so that, as shown in equation (51) 

H m ,o(j) = Qm.oiO) + Qm,oi(j + KM) + Q m ,oi(j + 2K/4) + Q m , 0 i(j + 3K/4) 
H m ,iO) = Qm.oi(j) + Q m ,Q2(j + K/4) - Q m , 01 (j + 2K/4) - Q m , 0 i(j + 3K/4) 

19 
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H m , 2 (j) = Qm,oi(j) " Q mt oiO + K/4) + Q m(0 iO + 2K/4) - Q m , 01 G + 3K/4) 
H m , 3 (j) = Q m .oi(j) - Qm.020 + K/4) - Q mi01 (j + 2K/4) + Q m , 02 (j ♦ 3K/4) 
H m ,oG + K/4) = Q mtO2 0) + Qm t o20 + K/4) + Q m(02 (j + 2K/4) + Q m , O2 0 + 3K/4) 
H m ,iO + K/4) = Q m(O2 0) - Q m ,oiG + K/4) - Q m , 02 G + 2K/4) + Q m ,oi(j ♦ 3K/4) 
5 H m , 2 G + K/4) = Q m ,oi(j) - Qm,o 2 G + 2K/4) + Q mi02 (j + 2K/4) - Q mt02 {j + 3K/4) 

H m , 3 G + K/4) = Q m , 02 G) + Qm,oiG + K/4) - Q m(02 G ♦ 2K/4) - Q m ,oiG + 3K/4) 
0<j<K/4-1 (51) 
and by using the equation (34), an equation (52) 

10 X/2-1 

15 0 < k < K/4-1 , p = 0, 1 , 2, 3 - • ■ (52) 

holds, so that the volume of the integration and summation operations of equation (40) may be reduced to a 
quarter. 

It is seen from above that, with the method of the present invention, both MDCT defined by the equation 

20 (19 and I MDCT defined by the equation (3) may be realized by integration of a k number of terms and sum- 
mation, so that outputting of an N/2 number of signals may be executed by 

N * N/(2K) + 6N + KN/4 
number of times of integration and summation operations and a 

K * K/2 + N + N * N/(8 * K * K) 

25 number of coefficient tables, so that the MDCT and I MDCT operations may be realized by a simpler hardware 
even although the number of times of the processing operations and the coefficient tables are increased as 
compared to the conventional method. 

Returning to Figs. 2 ff., a preferred embodiment for practicing the above-mentioned principle of the illu- 
strated embodiment of the present invention is explained. 

30 That is, a signal transforming device for MDCT (MDCT calculation circuit 30) of the illustrated embodiment 
comprises an x 01 calculating circuit 33, as a forward transforming window unit for multiplying an input signal 
with a forward transforming window, and circuit components inclusive of an Xq 2 calculating circuit 34 as a linear 
forward transforming unit for linear forward transforming an output signal of the x 01 calculating circuit 34, and 
is adapted for processing an N number of samples of input signals. The linear forward transforming unit includes 

35 an X02 calculating circuit 34 and an Xo 3 calculating circuit 35, as pre-processing units for pre-processing an out- 
put signal of the X01 calculating circuit 35, and an integration and summation circuit 40 for processing an output 
signal of the pre-processing unit with integration and summation operations. The integration and summation 
circuit 40 executes integration and summation on an N/2 number of the input signals supplied from the pre- 
processing unit. The integration is executed on each of a K number of terms of the input signals, grouped to- 

40 get her as a processing unit for integration. The integration and summation operations are iteratively performed 
an N/(2*K) number of times for outputting a sum total of N/2 signals. 

In a flow chart of Fig.4, the signal transforming sequence in the embodiment of Fig.2 for executing MDCT 
operations by the method of the present invention is shown. At a first step S1 in Fig.4, chronological sample 
data, such as PCM audio data, are grouped into blocks each made up of a predetermined number of, herein 

45 N samples. The blocks are set so that an overlap with neighboring blocks amounts to 50%, that is the neigh- 
boring blocks are overlapped by N/2 samples, as described in connection with Fig.1, and sample data of the 
respective blocks are multiplied by a forward transforming window Wh, as described in connection with Fig.1. 
These filtered N-samples of the chronological data are modified at the next step S2 as shown by the equa- 
tions (3) and (4) for transformation into an N/2 number of real number data x 03 . At the next step S3, calculations 

so shown by the equations (17) to (20), (29) and (31) are executed on X03 to find the N/2 number of spectral data 
Xq which are outputted as MDCT output data. 

The above-described transformation for MDCT is executed by a signal transforming device for MDCT of 
the present embodiment, an illustrative hardware for realizing the transformation for MDCT, shown in Fig.2. In 
the signal transforming device, shown in Fig.2, chronological sample data x 0 , stored in chronological sample 

55 buffer 32, is transformed by an MDCT calculating circuit 30 for executing the signal transforming method for 
MDCT into spectral data which is outputted at an output terminal 39. 

Specifically, the MDCT calculating circuit 30 comprises an x 01 calculating circuit 33 for executing a proc- 
essing on the chronological sample data Xo according to equation (2) to find the chronological data X01, an Xq 2 
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calculating circuit for executing calculations according to equation (3) on the chronological data Xoi to find data 
Xo 2 , an calculating circuit 35 for executing calculations according to equation (4) on the data x 02 to find data 
Xoa, and an integration and summation circuit 40 for executing calculations according to equations (17) to (20), 
equations (33) and (35) on the data x 03 to find spectral data x<). 

5 Specifically, the integration and summation circuit 40 comprises an A^oi calculating circuit 41 for executing 

a processing according to equation (1 7) using data xoa and an output of an X calculating circuit 46 as later ex- 
plained to find data A^oi, an A,^ calculating circuit 42 for executing a processing according to equation (18) 
using the data x 03 and the output of the X calculating circuit to find data A m0 2. a B m0 i calculating circuit 43 for 
executing a processing according to equation (19) to find data B mi01 , a B m)02 calculating circuit 44 for executing 

10 a processing according to equation (20) to find data B mf02 , a G m calculating circuit 45 for executing a processing 
according to equation (33) using the data B mi01 and data B m 02 to find data G m , and an x 0 calculating circuit 46 
for executing a processing according to equation (35) to find spectral data xq. 

Meanwhile, the xqi calculating circuit 33 as a forward transforming window unit, and the xq 2 calculating 
circuit 34, as a pre-processing unit for the linear forward transforming unit, may be grouped into one circuit. 

15 A signal transforming device for IMDCT according to an embodiment of the present invention (IMDCT cal- 

culating circuit 50) comprises a linear inverse transforming unit for executing linear inverse transform on an 
input signal and an inverse transforming window unit for multiplying an output signal of the linear inverse trans- 
forming unit with an inverse transforming window, as shown in Fig.3, and is adapted for performing inverse 
transform operations on an N/2 number of independent input signals. The nonlinear inverse transforming unit 

20 comprises a y u calculating circuit 52 for pre-processing an input signal, an integration and summation circuit 
60 for executing integration and summation on an output signal of the y n calculating circuit 52, and a Wi cal- 
culating circuit 56 for post- processing an output signal of the integration and summation circuit 60. The inte- 
gration and summation circuit 60 executes integration and summation operations on the N/2 number of input 
signals supplied from the y u calculating circuit 52 by grouping the input signal into groups each made up of K 

25 terms, and iteratively performs the operations N/(2*K) times to output a sum total o N/2 number of signals. 
Meanwhile, in the example of Fig.3, the Vv^ calculating circuit 56 is composed of the inverse transforming unit 
and a post-processing unit of the linear inverse transforming unit. 

In a flow chart of Fig.5, the signal transforming sequence in the embodiment of Fig.3 for executing IMDCT 
operations by the method of the present invention is shown. At a first step S11 shown in Fig.5, an N/2 number 

30 of spectral data y 1 is processed by sign exchange and re-arraying as shown by equation (9) for being converted 
into an N/2 number of real number data y^. At the next step S 12, the real number data y^ are processed by 
processing operations in accordance with equations (32) to (35), (43) and (44) to generate an N/2 number of 
data W 0 . At the next step S15, the data W 0 are processed by transform processing as shown by equation (37) 
to find an N number of real number data Wt. 

35 The above-described transform operations for IMDCT and the transform processing for constructing 

chronological data from the output of the IMDCT transform operations are executed by a signal transforming 
device for MDCTof the present embodiment, an illustrative hardware for realizing the transformation for IMDCT, 
shown in Fig.3. The spectral data y-i from input terminal 51 are multiplied by an inverse transform window by 
a chronological data buffer 57 and subsequently summed to IMDCT output data of both neighboring blocks by 

40 an overlap addition circuit 58 before being outputted as chronological sample data at output terminal 59. 

Specifically, the IMDCT calculating circuit 50 comprises a y n calculating circuit 52 for executing calcula- 
tions on the spectral data y^ in accordance with the equation (9) to find data y 11p an integration and summation 
circuit 60 for executing calculations on data y t1 in accordance with the equations (36) to (39), (51) and (52) to 
find data W 0 , and a W t calculating circuit 56 for executing calculations on data W 0 in accordance with the equa- 

45 tion(41)tofinddataW v 

On the other hand, the integration and summation circuit 60 comprises a P mi01 calculating circuit 61 for 
executing calculations on the data y„ and an output of a W 0 calculating circuit 66 as later described to find 
data P mt0 i, a P mi02 calculating circuit 61 for executing calculations on the data y^ and the output of W 0 calcu- 
lating circuit 66 as later described to find data P mr0 2> a Qm,oi calculating circuit 64 for executing calculations on 

so the data P mt0 i and P mr02 in accordance with the equation (38) to find data Q m , 0 t. a Q m ,o2 calculating circuit 64 
for executing calculations on the data P m0 i and P m0 2 accordance with the equation (39) to find data Q mt0 2» 
an H m calculating circuit 65 for executing calculations on data Q mi0 i and Q m>0 2 in accordance with equation (51) 
to find data H m and a W 0 calculating circuit 66 for executing calculations on the data H m in accordance with 
equation (52) to find data W 0 . 

55 Referring to Fig.6, explanation is made of a concrete example of a high efficiency encoding device, to which 

MDCT as one of the above-described signal transforming operations of the illustrated embodiments is applied. 
The concrete high efficiency encoding device shown in Fig.6 makes use of techniques of sub-band coding, 
adaptive transform coding and adaptive bit allocation, That is, with the high efficiency encoding device, shown 



21 



EP 0 590 790 A2 



in Fig. 6, input digital signals, such as input PCM audio signals, are split into plural frequency bands so that the 
bandwidth becomes broader towards the higher frequency. The band-limited signals are processed by MDCT, 
a kind of orthogonal transform processing operations, for each frequency band, to generate frequency-domain 
spectral data, which are then encoded by adaptive bit allocation for each of so-called critical bands. 

5 That is, referring to Fig.6, audio PCM signals in e.g. a range of from 0 to 20 kHz are supplied to input ter- 

minal 11. The input signals are split by a band-splitting filter 12, such as a so-called QMF filter, into a band of 
0 to 10 kHz ad a band of 10 to 20 kHz. the band of 0 to 10 kHz is split by a similar band-splitting filter 13, such 
as QMF filter, into a band of 0 to 5 kHz and a band of 5 to 10 kHz. The signals in the band of from 10 kHz to 
20 kHz from band-splitting filter 13 are transmitted to the modified DCT (MDCT) circuit 14 of the illustrated 

10 embodiment for MDCT processing operations. Similarly, the signals in the band of 5 to 1 0 kHz and in the band 
of 0 to 5 kHz, supplied from the band-splitting filter 1 3, are transmitted to MDCT circuits 15, 16, respectively, 
for MDCT processing operations. 

The spectral data on the frequency domain or coefficient data, generated by MDCT processing operations 
by the MDCT circuits 14 to 16, are grouped for each of the critical bands so as to be supplied to an adaptive 

15 bit allocation and encoding circuit 17. 

Meanwhile, the critical band means a frequency band splitting system which takes human acoustic sense 
into consideration. That is, the critical band is a band of a narrow bandwidth noise which masks a pure tone 
having a frequency in the vicinity of that of the noise. The critical band becomes broader in bandwidth towards 
a higher frequency. The entire audible frequency range of from 0 to 20 kHz is split into e.g. 25 critical bands. 

20 The adaptive bit allocation and encoding circuit 1 7 normalizes each of the spectral signals by a scale factor, 

that is a maximum value of the absolute values of the spectral signals contained in each critical band, from 
one critical band to another. The circuit 17 then re-quantizes the normalized spectral signals with a sufficient 
number of bits to mask the quantization noise by the critical band signals to output the re-quantized spectral 
signals along with the scale factor as found for each critical band and the number of bits employed for requan- 

25 tization. The data encoded in this manner are outputted via output terminal 18. 

Referring to Fig. 7, explanation is made of a concrete example of a high efficiency decoding device, to which 
IMDCT as one of the above-described signal transforming operations of the illustrated embodiments is applied. 

Referring to Fig. 7, the above-mentioned scale factor, the number of bits employed for requantization and 
requantized and encoded spectral signals are entered at input terminal 21 and spectral signals are generated 

30 by a spectral signal decoding circuit 22 from these input data. Of these spectral signals, the spectral signal in 
the ranges of 10 to 20 kHz, 5 to 10 kHz and 0 to 10 kHz are transformed into signal waveform data of the 
respective bands by IMDCT circuits 23, 24 and 25, respectively. Of these signal waveform data in the three 
frequency bands, the signal waveform data in the band of 0 to 5 kHz are combined with the signal waveform 
data in the band of 5 to 10 kHz by a band combining circuit 26 into signal waveform data in the band of 0 to 

35 10 kHz which is then combined with the signal waveform data in the band of 20 to 20 kHz by a band combining 
circuit 27 into signal waveform data over the entire frequency range which is outputted at output terminal 28. 

It should be noted that the present invention is not limited to the above-described embodiments. For ex- 
ample, the device of the present invention may be appl ied to a variety of transform encoding circuit or a decoder 
in addition to the high efficiency encoding/decoding device shown in Figs.6 and 7. 

40 The forward transform calculations for modified DCT (MDCT) may be performed efficiently on an n number 

of chronological sample data sliced with 50% overlap with both neighboring blocks, using FFT with a length 
of N/4, for reducing the processing volume or work area for constructing a high efficiency encoding device for 
audio signals. Besides, inverse transform calculations for modified DCT (MDCT) may be performed efficiently 
on an N/2 number of independent spectral data using FFT with a length of N/4 for reducing the processing 

45 volume or work area for constructing a high efficiency decoding device for audio signals. 



Claims 

50 1. A modified DCT signal transforming device for processing an N number of sample input signals, compris- 
ing a forward transforming window unit for multiplying input signals with a forward transforming window 
and a linear forward transforming unit for linear forward transforming output signals of said forward trans- 
forming window unit, characterized in that said linear forward transforming unit comprises a pre-process- 
ing section for pre-processing the output signals of said forward transforming window unit and an inte- 

55 gration and summation processing section for executing integrating and summation processing operations 

on output signals of said pre-processing section, said integration and summation processing section exe- 
cuting integration and summation processing operations on an N/2 number of input signals supplied from 
said pre-processing section, said N/2 number of the input signals being previously grouped into sets each 
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consisting of a K number of terms, said integration and summation operations being iteratively executed 
an N/(2*K) number of times for outputting a sum total of N/2 number of signals. 

The modified DCT signal transforming device as claimed in claim 1 wherein K in said integration and sum- 
mation processing section is a power of 2 which is not less than 4. 

The modified DCT signal transforming device as claimed in claim 1 wherein, with an N number of samples 
of input signals of from 0 to (N-1), obtained on multiplying the forward transforming windows in the forward 
transforming window section, the pre-processing unit outputs the n+3N/4th input signal inverted in sign 
as the nth intermediate signal for n of from 0 to N/4-1, and an n-N/4th input signal as the nth intermediate 
signal for n of from N/4 to N-1, and wherein, with an N number of samples of intermediate signals of from 
0 to (N-1), the pre-processing unit outputs the 2n-th intermediate signal less the N-1-2n-th intermediate 
signal as an nth output signal. 

The modified DCT signal transforming device as claimed in claim 1 wherein the integration and summation 
processing unit decreases the number of times of the processing operations by grouping an N/2 number 
of input signals supplied from said pre-processing unit into sets each consisting of a K number of terms, 
based on periodicity, for outputting said N/2 number of the signals. 

The modified DCT signal transforming device as claimed in claim 1 wherein said integration and summa- 
tion processing unit performs, on an N/2 number of input signals supplied from said pre-processing unit, 
a series of processing operations for m of from 0 to N/(2K)-1 until the number of output signals becomes 
equal to N/2, said series of the processing operations consisting in providing, as a jth first intermediate 
signal, a sum for Is of from 0 to N/2-1 of a product of the (Kl+j)th input signal multiplied by 
cos(rc(2m+1)(4K+1)/(2N), for j of from 0 to k-1, providing, as a jth first intermediate signal, a sum for Is of 
from 0 to N/2-1 of a product of the (K(l-1)+j)th input signal multiplied by sin(rc(2m+1)(4KI+1)/(2N), for j of 
from K to (2K-1), providing, as a jth second intermediate signal, a sum of the jth first intermediate signal 
multiplied by -cos(2rc(2m+1)(j/N) and the (j-K)th first intermediate signal multiplied by - sin(2n(2m+1 )(j/N), 
for j of from K to K-1, providing, as a jth second intermediate signal, a sum of the (j-K)th first intermediate 
signal multiplied by -sin(2n(2m+1)(j-K)/N) and the jth first intermediate signal multiplied by -cos(2n(2m+1)(j- 
K)/N), for j of from K to 2K-1, providing, as a k-th third intermediate signal, a sum of the jth second inter- 
mediate signal multiplied by cos(2nJk(4j+1)/(4K)), for j of from 0 to K-1, providing, as the k-th third inter- 
mediate signal, a sum of the jth second intermediate signal multiplied by sin(2nk(4(j-K)+1)/(4K)), for k of 
from K to (2K-1), and providing, as a (m+kN/(2K))th output signal, a sum of the k-th third intermediate 
signal and a (m+Kn/(2K))th third intermediate signal, for k of from 0 to K-1. 

The modified DCT signal transforming device as claimed in claim 5 wherein said integration and summa- 
tion processing unit performs a processing operation of computing an output signal of the integration and 
summation processing unit from the second intermediate signal by a processing volume reduced by one 
half by grouping the second intermediate signals by each two terms by taking advantage of periodicity 
so that the number of output signals becomes equal to N/2. 

The modified DCT signal transforming device as claimed in claim 5 wherein said, integration and sum- 
mation processing unit performs a processing operation of computing an output signal of the integration 
and summation processing unit from the second intermediate signal by a processing volume reduced by 
one quarter by grouping the second intermediate signals by each four terms by taking advantage of per- 
iodicity so that the number of output signals becomes equal to N/2. 

A modified DCT signal processing device for processing an N/2 number of independent input signals, com- 
prising a linear inverse transforming unit for linear inverse transform of input signals, and an inverse trans- 
forming window unit for multiplying an output signal of said linear inverse transform unit by an inverse 
transforming window, 

said linear inverse transforming unit comprising a pre-processing unit for pre-processing an input 
signal, an integration and summation processing unit for executing an integration and summation proc- 
essing operations on an output signal of said pre-processing unit, and a post-processing unit for post- 
processing an output signal of said integration and summation processing unit, said integration and sum- 
mation processing unit executing the integration and summation operations on the N/2 number of input 
signals supplied from said pre-processing unit by grouping the input signal by a K number of terms, and 
executing the operations iteratively by N/(2*K) number of times for outputting a sum total of an N/2 number 
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of signals. 

9. The modified DOT signal transforming device as claimed in claim 8 wherein K in said integration and sum- 
mation processing section is a power of 2 which is not less than 4. 

10. The modified DCT signal transforming device as claimed in claim 8 wherein said pre-processing unit out- 
puts a signal equivalent to a k-th output signal which is the 2k-th input signal for values of k of from 0 to 
N/4-1 and a signal equivalent to a k-th output signal which is the (N-1-2k)th input signal inverted in sign 
for values of k of from N/4 to N/2-1 . 

11. The modified DCT signal transforming device as claimed in claim 8 wherein said integration and summa- 
tion processing unit decreases the number of times of processing operations for integration and summa- 
tion processing of N/2 number of input signals supplied from said pre-processing unit based on periodicity 
for outputting said N/2 number of signals. 

12. The modified DCT signal transforming device as claimed in claim 8 wherein said integration and summa- 
tion processing unit performs, on an N/2 number of input signals supplied from said pre-processing unit, 
a series of processing operations for m of from 0 to N/(2K)-1 until the number of output signals becomes 
equal to N/2, said series of the processing operations consisting in providing, as a jth first intermediate 
signal, a sum for Is of from 0 to N/2-1 of a product of the (Kl+j)th input signal multiplied by 
cos(7i(2m+1)(4K+1)/(2N), for j of from 0 to K-1, providing, as a jth first intermediate signal, a sum for Is of 
from 0 to N/2-1 of a product of the (K(l-1)+j)th input signal multiplied by sin(*(2m+1)(4KI+1)/(2N), for j of 
from K to (2K-1), providing, as a jth second intermediate signal, a sum of the jth first intermediate signal 
multiplied by-cos(2n(2m+1)(j/N) and the (j-K)th first intermediate signal multiplied by - sin(2n(2m+1)(j/N) t 
for j of from K to K-1, providing, as a jth second intermediate signal, a sum of the (j-K)th first intermediate 
signal multiplied by -sin(2rc(2m+1)(j-K)/N) and the jth first intermediate signal multiplied by -cos(2n(2m+1)(j- 
K)/N), for j of from K to 2K-1 , providing, as a k-th third intermediate signal, a sum of the jth second inter- 
mediate signal multiplied by cos(2nJk(4j+1)/(4K)), for j of from 0 to K-1, providing, as the k-th third inter- 
mediate signal, a sum of the jth second intermediate signal multiplied by sin(27tk(4(j-K)+1)/(4K)), for k of 
from K to (2K-1), and providing, as a (m+kN/(2K))th output signal, a sum of the k-th third intermediate 
signal and a (m+kN/(2K))th third intermediate signal, for k of from 0 to K-1. 

13. The modified DCT signal transforming device as claimed in claim 12 wherein said integration and sum- 
mation processing unit performs a processing operation of computing an output signal of the integration 
and summation processing unit from the second intermediate signal by a processing volume reduced by 
one half by grouping the second intermediate signals by each two terms by taking advantage of periodicity 
so that the number of output signals becomes equal to N/2. 

14. The modified DCT signal transforming device as claimed in claim 12 wherein said integration and sum- 
mation processing unit performs a processing operation of computing an output signal of the integration 
and summation processing unit from the second intermediate signal by a processing volume reduced by 
one quarter by grouping the second intermediate signals by each four terms by taking advantage of per- 
iodicity so that the number of output signals becomes equal to N/2. 

15. The modified DCT signal transforming device as claimed in claim 8 wherein said pre-processing unit out- 
puts a signal equivalent to an (n+N/4)th input signal for the value of n of from 0 to N/4-1, a signal equivalent 
to a 3N/4-1-nth input signal reversed in sign for values of n of from N/4 to 3N/4-1 and a signal equivalent 
to an n-3N/4th input signal reversed in sign for values of n of from 3N/4-1 to N/4-1 . 
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